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Abstract: We discuss the general decomposition and possible general parameterizations 
of the processes e^e~ ^ 7* — > P1P27, where P1P2 = tt^tt^ , vr'^r/, or vr+vr^, for ^/s ^ M$. 
Particular attention is devoted to the amplitude where the two pseudoscalar mesons are 
in a J*^^ = O^"*" state, where we propose a general parameterization which should help to 
shed light on the nature of light scalar mesons. 
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1. Introduction 

As widely discussed in the literature, radiative $ decays, such as ^ /o7 tttt'j or 
^ Oo7 ^ '^Vl: of the primary sources of information about the interesting and 

still controversial sector of light scalar mesons (see e.g. Refs. ^, ^ and references therein). 
In principle, the large amount of data collected at <I> factories should allow to study these 
processes with excellent accuracy. However, it must be realized that these processes are 
only one of the components of the basic ^-factory observables, namely the e~^e~ — > -P1P27 
cross sections (where P1P2 = n^ir^, vr'^r/, or 7r"^7r~). An accurate and possibly model- 
independent description of all the components of these reactions is necessary in order to 
extract reliable information about the scalar sector of QCD. 

The purpose of the present paper is twofold. First, we discuss how to isolate the 
interesting scalar amplitude from the other contributions to the cross section. Second, we 
present a general parameterization of the scalar amplitude which should allow a model- 
independent determination of basic parameters such as masses, widths and couplings of /o 
and ao mesons. 

The paper is organized as follows. In Section |2| we discuss the general decomposition 
the cross sections in terms of gauge-invariant tensors. We start from the simpler case 
of neutral final states, |7r''7r^7) and |7r^?77). We then generalize to the |7r+7r^7) case, 
where initial- and final-state radiation represent a serious background. In Section |^ we 
analyse the contributions to the invariant amplitudes due to the exchange of vector mesons 
(e+e" — > VPi(^2) ~^ P1P27), which represent an irreducible physical background for the 
scalar amplitude. Finally, in Section ^ we present a general decomposition of the scalar 
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amplitude (e^e^ — > 5*7 — > P1P27) which takes into account the narrow-width structure 
of /o and oq at high Mp^p^, and the constrains of unitarity and chiral symmetry at low 
Mp^Pj. The complete procedure we propose for the fit of the cross sections is summarized 
in the last section. 



2. General decomposition of the cross sections 
2.1 Neutral final states 

The generic matrix element for the transition e+e^ Pj'PgT (-^1,2 = or r]) can be 
written as 



where s = (p+ = P^- The constraints of gauge invariance imply 



P^T^"' = KT^" = . 



(2.1) 



(2.2) 



Using this notation, the differential cross section with unpolarized beams becomes 



spins 



— nP^ T T* 



where 



Spi d^p2 



.(27r)*5^^\P-pi-p2-k) 



, (2.3) 



(2.4) 



(27r)32£'^ (27r)32£;i (27r)32^2 
and C12 takes into account the 1/2 factor in case of identical particles: C12 = 1 for 
{P^PO} = {V^h Cu = 1/2 for {P0,P0} = {^o,vrO}. 

As a result of the gauge-invariance conditions, in the case of a real photon in the final 
state we can decompose T^'^ as the sum of three independent structures: 



T, 



AiL«+^24'J+A34'2 



(2.5) 



with 



L^} = ^i(P- ^) [(^ • Q)9,.u - k^qA + {k ■ P)qf,qu - {k ■ q)qf,Pu} 
= -[ik- q){P^gt.u - P^Pu) + {k ■ P)Pf,qu - P'k^q,] , 



q = ^(p^ -P^) 



(i) 

The normalization of T^jy and the -L^,> has been chosen in order to maximize the contact 
with Ref. 0: Lji^ and Lg] are identical to the corresponding expressions of Ref. while 
L^fj] coincides with the corresponding tensor of Ref. |Q] only in the case of identical particles. 
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The dimensionless form factors Ai are determined by the specific dynamical model. 
For instance, the scalar amplitude e~^e~ — > S'j — > -P1-P27 induces a non- vanishing contri- 



bution only to Ai, while the transition e~^e~ — > VP^ 
contributions to all the form factors. 

The phase space element can be re-written as 



1(2) 



P1P27 leads to non- vanishing 



1 



-dEiclE^dni 



1 



.dEidE2dnb , 



where -£^7,1,2 are the energies in the c.o.m. frame and fib denotes the solid angle of the 
beam axis with respect to the decay plane. ^ Performing the angular integration we find 



F{x,xi,X2) = ^ / dr^b 



— nP'^T T* 



aii|Ai|2 + a22|^2p + assl^sl^ + ai2{AiAl + A\A2) 
+ai^{AiAl + A\Az) + 023(^2^^ + AlAz) , 



where the aij coefficients are given by 



an = 


^ 2 

-X 




3 






ai2 = 


|[(.. 






ai3 = 


-x{xi 






0^22 = 


>^ 




-5{xi 



0-23 
^33 



{{xi - X2f + 2x^{a -l + xf + 2(xi - X2f{l - x){a - I + x) 

5{xi - X2) [2{xi - X2f + 2(0- - 1 + x){xi + X2)] 
[{x^-x2f + 2{a -l + x)]] 

^(xi - X2) [{xi - X2f - 25{xi - X2) + 5^] 

- [(xi - X2)^(l + x)- x^{a - 1 + x) - 5{xi - X2)(x + 2) + J^] 



in terms of the adimensional variables 

2Si 

rs 

r2 I ^j2\ 



2E^ 



Xi 



a 



2{Mf + M|; 



2 - Xi - X2 

2(Mf - M|) 



coincide with the Cjj of Ref. 



s s 

In the limit 6^0 the aij coefficients in Eq. (| 

the exception of ai2, which is reported incorrectly in Ref. [Q]. 

The general expression for the total cross-section is then given by 



a 



327r2s 
a 

32tt^s 



C12 F{2 — xi — X2, xi, X2) dxidx2 
C12 F{x, xi, 2 — X — xi) dx dxi. 



(2.7) 



(2.8) 



(2.9) 
with 



(2.10) 



^ The element of the sohd angle can be expressed as df2b = dcos6jd(f>, where 6j is the angle between 
photon and momenta, and (/> the orthogonal angle with respect to the decay plane. The latter leads to 
a trivial integration in the \Ai\^ term, but is non-trivial for the other contributions. 
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c) [fsr; 



Figure 1: Tree- level diagrams for ISR, FSR and SD contributions to e+e 
amplitude corresponds to the 0{E^) gauge-invariant terms in the diagram e). 



TT+TT 7. 



The SD 



2.2 Initial- and final-state radiation and generalization to the e+e — > vr+vr 7 
case 



TT^TT 7: 



At 0{a^) we can distinguish three basic components to the processes e+e" 
initial-state radiation (ISR), final-state radiation (FSR) and the irreducible structure- 
dependent (SD) amplitude which vanishes in the ^ Q limit (see Fig. ||). The latter, 
which is identical to the neutral case discussed before, is the only contribution sensitive to 
scalar-meson dynamics. 

By construction, ISR and FSR amplitudes can be fully described in terms of the elec- 
tromagnetic form factor of the pion and are known with good accuracy^ (for an extensive 
discussion applied to the ^-factory case see e.g. Ref. ||6|, 0]). The Ivr^vr") state produced 
by ISR has opposite transformation properties under parity and charge conjugation with 
respect to those produced by FSR and SD. For this reason, as long as the kinematical cuts 
applied on the cross section are symmetric under the exchange vr"*" <-> vr", we can neglect 
the interference of the ISR amplitude with the other two contributions. On the other hand, 
the FSR amplitude, which can be decomposed according to Eq. ( |2.5| ), does interfere with 
the SD terms: 



d(T(e+e vr+TT 7) = duisR + dcrpsR + d^-SD + f^f^rsR+SD • 



(2.11) 



The form factors of the FSR contribution, which are singular in the x ^ limit, can be 
written as 



^FSR 



47raj;F,r7r('S) 
— (xi — X2Y 



/.FSR 
^2 



— {xi — X2Y ' 



Af"^ = , 



(2.12) 



^ We denote here as FSR amplitude only the gauge-invariant part of the amplitude which can unam- 
biguously related to the on-shell non-radiative e+e~ — > ■k'^-k~ processes, according to Low's theorem |5|. 
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Figure 2: Leading contributions to the structure-dependent amplitude: vector-meson exchange 
(left) and scalar (right) contributions. 

where -F7r7r(s) denotes the pion electromagnetic form factor (with the standard normaliza- 



tion F7r7r(0) = 1). The last three terms in Eq. ( 2.11 ) can thus be described by a generic 



expression of the type in Eq. ( 2.10| ), where the form factors are obtained by summing the 



singular ^^'^^ to the regular terms of the SD amplitude. 

So far we have analysed only the ©(a^) contributions to the cross sections, or the lead- 
ing contributions in the case of a single detected photon in the final sate. As is well know, 
the physical cross sections are obtained by the convolution of these leading expressions with 
appropriate radiation functions |^ which take into account the effect of the undetected ISR 
(both for charged and neutral final states) and the undetected FSR (for the charged final 
state only). 

3. The vector-meson exchange contribution 

In a good fraction of the allowed kinematical region, a sizable contribution to the structure- 
dependent amplitude is induced by the vector- meson exchange process in Fig. § (left). 

The starting point to evaluate this contribution are the effective couplings 5^^, and 
ffp^ (with dimension 1/energy), defined by 

M [V{e,P) ^ P{qh{e,k)] = egl^ e'^^f'^e^e.q.K , 
M [V{e, P) ^ P{q)V'{e, k)] = g^^, e^-'^'^e^e.q.k^ , (3.1) 

and the adimensional vector- meson electromagnetic couplings, Fy, defined by 

M [V{e) ^ e+e-] = ^ u{p^)^^v{p+) . (3.2) 

In terms of these couplings, the vector-meson exchange process in Fig. ^ give rise to the 
following contributions to the form factors: 

AT' = -^-1 + + ^) ^9{xi) + g{x2)\ + \[xi-X2 + [g{xi) - g{x2)] 

Al^^ = \[g{x,)+g{x2)] 

AT' = -\[9{x,)-g{x2)] (3.3) 



where 



- 5 - 





i'P7 L J 




\Fp\ = 16.6 ±0.2 


Ig" 1 = 0.79 ±0.10 

— i.OD It u.uo 


|o" 1 = 12 ± 2 

1 1 1 o _|_ n 9 

u — i.o zt U.Z 

= (4.1 ±0.5) X 10~2 
|o* 1 = 1.1 ±0.2 

In* 1 = 210" 1 = 21 0" 1 = 7.0 ± 1.5 


= 55.9 ±0.5 


|o- 1 = 2.38 ± 0.04 
\gp 1 = 0.46 ±0.02 


|Fci,| = 44.4 ±0.4 


|o* 1 = 0.132 ±0.004 
|5* 1 = 0.691 ± 0.008 



Table 1: Reference values for the modulos of the effective vector-meson couphngs defined in 



Eqs. (ly)-(^ 



and 



Dx{q^)=s-Ml+iMxTx 



(3.5) 



In the limit of identical particles in the final state, these results are fully consistent with 
those of Ref . Q . As discussed in Ref . , a natural improvement of the above expressions 
is obtained by the replacement of the constant Tx with appropriate energy-dependent 
widths, taking into account the velocity factors of the dominant final states 

In Table |l] we report the current estimates for the most relevant set of Fy, g"^ , and 



9' 



^, couplings. The results for Fy and g^_^ have been determined by means of the relations 



r(y 



e ' e 



aMv 



V{V ^ P7) 



a\g 



P7I 



M2 



M 



2n3 



2Mv 



(3.6) 



using the experimental values of r(y e^e^) and T{V — > P7) in |10|. The 5^^, have been 
determined theoretically,^ with the exception of 5* 
from the experimental value of r(<I> 



which has been determined directly 
wvr) in |1C]. We stress that the results reported in 
Table |l| should be considered only as rough reference values, or as natural starting point 
for a fit of the cross sections. The high-statistics data on the e^e" P1P2I reactions at 
a $ factory should allow to determine these couplings (or at least some combinations of 
them) with much higher accuracy. 



4. The scalar amplitude 



We are now ready to analyse the scalar amplitude or, more precisely, the contributions to 
the form factor Ai not described by vector mesons and/or FSR. 

The contribution to Ai induced by the exchange of a single vector and a single scalar 
resonance, as shown in Fig. § (right), is: 



^i(e^ 



V ^S-i^ P1P27) 



4Fy Dv{s)Ds\{l-x)s\ 



(4.1) 



^ The g^^i have been determined by: i) assuming a simple effective Lagrangian of the type 
C = g tr[{V,V'}P), where V, V' , and P are 3x3 matrices in flavour space; ii) fixing the $-0; mix- 
ing angle from the Fv values; iii) enforcing the vector-meson dominance relation = X^y' 9^^' I 
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where the couplings gf^ and are defined by 



T{S ^ P1P2 



a\g 



S7I 



M2 -M| 



2My 



(4.2) 
(4.3) 



[M2 - (Ml - M2)2] [M2 - (Ml + M2) 
2M 



211/2 



If we could consider only this resonant contribution, the total cross-section would assume 
the following form 



das. 



icalar 



^Ci2 \Ai{e+e- ^V^S-f^ PiP2-/)\^ E^dE^dEi 



BWv{s[ 



r{V e+e') 
MvTl 



12^57 



^5(^12) 



dE, 



'S12 



7 ' 



where 



(4.4) 



Sl2 



s — 2\[sE^ = (1 — x)s 



\Dv{s)\^ (M 



(4.5) 



Note that the Et factor in (4^), which is dictated by gauge- invariance according to the 



general decomposition (2.8), implies a sizable distortion from a standard Breit-Wigner 
shape for resonances close to the end of the phase space, such as /o(980) and ao(980). 

The simplified expression ( |4.4| ) is often used in the literature to describe the contribu- 
tions of the narrow resonances /o(980) and ao(980). However, a coherent description of all 
the amplitudes contributing to the physical processes requires a more refined treatment. 
First, in order to compute the total cross section we need to consider the general expression 
in Eq. ( 2.10| ), with coherent sum of all contributions to the Ai'. 



full 



^FSR _j_ ^vcct _|_ 



seal 



yj^FSR _|_ ^vect 



A 



full 



/I vect 
^3 



(4.6) 



with Af^^ (tt+tt" case only) and jd™'^* given in Eq. ( p. 121 ) and Eq. p.3| ), respectively. 
Second, the expression of J^^^ can involve several resonances and, possibly, also non- 
resonant backgrounds. 

Since /o(980) and ao(980) are the only two narrow scalar resonances with mass below 
M<j>, a convenient parameterization for J^^^ is 



^scal 



E 



sMl 



V=p,LO,(t>; S=fo,ao 



AFv Dv{s) 



Ds{s 



12 



+ <2(512) 



(4.7) 



Here i?i^(si2) denotes a non-resonant term which, in absence of re-scattering effects, can 
be expanded as a regular series in powers of (si2 — M|). The situation becomes particularly 
simple in the exact isospin limit, where a single narrow resonance can contribute to each 
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channel: the fo{I = 0) in the two Ivrvr) cases and the ao(/ = 1) in the |vrr/) one. In this 
limit, at fixed values of s and neglecting re-scattering effects, we expect a structure of the 
type 

with a different set of effective couplings for each channel. As we will discuss in the 
following, this structure can be systematically improved in order to take into account the 
elastic re-scattering phases of the two pseudoscalar mesons. 

It is worth to stress that a sizable non-resonant term is phenomenologically required 
from data, at least in the |vrvr) channels. Indeed, if we only retain the pole term in Eq. ([4.^), 
the factor in the cross section -see Eq. ( |4.4| )- implies a too large result at low Stttt 
compared to observations (see e.g. Ref.|jll|, |l^): experimental data clearly indicate that at 
low Stttt the contribution of the /o(980) is partially compensated by other contributions. 

In the |vrvr) channels one can consider an alternative parameterization where, in addi- 
tion to the narrow /o(980), also the broad /o(600) is included by means of an appropriate 
complex propagator in the / = channel.^ In this case the /o(600) could be responsible for 
the partial cancellation of the /o(980) contribution at low Stttt, with a minor role played by 
the non-resonant term. In principle, precise enough experimental data on the e^e~ — > vrvr7 
cross sections should be able to distinguish the case of an explicit pole structure for the 
/o(600) from a pure polynomial term. However, the broad nature of the /o(600) makes 
this distinction very difficult, even with the high statistics available at a <I> factory. 



4.1 Re-scattering phases for |vrvr) final states 

As anticipated, the parameterization (|4.8| ) can be improved in order to take into account the 
absorptive parts generated by elastic re-scattering. We illustrate here how this improvement 
can be implemented in the |vrvr) case, where this effect is more relevant and the information 
about elastic re-scattering at low St^t^ is very precise. For simplicity, we consider only 
mediated contributions and we include only the /o(980) as explicit pole structure. 

Neglecting non-<I> contributions, we can decompose the scalar form factors for charged 
and neutral |vrvr) final states as 



Ai yn vr 



Af^l(vr°vr°) 



sMl 



77scal 
^2 



(St. 



4F$ D$(s) 
e sM| 
4F$ Dq,{s) 



%/3 



\/2F|-i(s.-) 



TTiscal 



(■^tttt) 



V2> 



(4.9) 



* The existence of a pole in the S-wave, 7 = 0, tttt scattering amphtude -corresponding to the /o(600)- 
is not under doubt However, this pole is very far from the real axis and quite close to the tttt 



threshold [g^. As a consequence, its contribution to the amphtude (4.7) is not necessarily well described 
by a simple complex propagator as in the /o(980) case. 
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where St^tt = (1 ~ x)s and the reduced F^^f^ correspond the / = 0, 2 isospin combinations. 
If we include only the /o(980) as explicit pole structure, we then have^ 



-iscal 



(^tttt) 



gfo gf 



Fl 



iscal / 



2 \S^^) = R2is^n) ■ (4.10) 

For Stttt close to the /o(980) pole, the leading elastic and inelastic re-scattering effects 
are automatically included in the scalar resonance propagator Q. In particular, considering 
only two-body intermediate states (tttt and KK) and defining 

S/o(s) = iV^rj.is) and Vi{s) = (s/4 - M^)'/^ , (4.11) 
the effective energy-dependent width assumes the form 
4M2 < s < 4M|± ^f,is) = ^ {-i9{^Kf[\vK+{s)\ + Ms)\\ + i{gll)^v^is)} , 



Svr^ 
i 



{(5xV)' l^K+i''^) + ^i^o(s)] + (gllfv^is)} . (4.12) 



The key observation which allows to determine the absorptive parts of the form factors 
in the low s-jrir region is the Fermi- Watson theorem jl^ . This implies that below the inelastic 
threshold the phases of the two Ff^^^lsT^T^) coincide with the vrvr S'-wave elastic phases. 
The low-energy structure of tttt phase shifts 
is known very precisely [15|. In the 5-wave 
channels we can write 



tan[5/(s)] 



ai + bi- 



So 





ai 


bi 


1 = 


0.220 ±0.005 


0.275 ± 0.009 


1 = 2 


0.045 ±0.001 


0.081 ± 0.002 



Table 2: Numerical values for the low-energy 
parameters of S'-wave tttt phases from Ref. |15| . 



with So = 4M^, and the aj and bj reported 
in Table |. 

In the 1 = 2 case, where the resonance term is absent, the Fermi- Watson constraint 
applies directly to the non-resonant term R2{s-k-k)- Generalizing the polynomial expansion 
in Eq. ([4.8|), this implies 



-^2(57. 



Ml) + 0[{s, 



„i(52(s,r7r) 



(4.13) 



where the 7^ are real parameters, as implicitly assumed for all the effective couplings so far 
introduced (g^^, 5^^, 5^ . . . ). 

The condition to be imposed on the / = term is slightly more complicated since 
we cannot ignore the phase shift induced by the /o propagator. Given that below the 
inelastic threshold all phase shifts are proportional to the pion velocity W7r(s,r7r), a convenient 
decomposition for i?o(svr7r) is 



-Ro(s7r7r) - 



(s„ 



Ml) + 0[[s,, 



(4.14) 



^ The effective coupling gi% denotes tiie coupling of the /o(980) to the two-pion state with 1 = 0. 
Correspondingly, the decay width appearing in the definition (4.3) must be interpreted as r[/o ~* (7r7r)/=o] = 

r(/o ^ 7 ■ 



-) + r(/o ^ttV"). 
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where at and /3j are real parameters. Imposing the Fermi- Watson constraint allows to 
determine the Pi in terms of the and the parameters of (5o(s^7r). Indeed, expanding 



the phase of the full form factor in powers of {s-„ 
following two conditions: 

Ml 



So), up to second order, leads to the 



(3o = ^[aso)-{so-Ml)C'{so)] 



ai 



(4.15) 



where 



+ 



{s - Ml + Re[S^„(s)])2 + Im[S/„(s)]2 ^ 



So 



So 



,_M|^+Re[S;„(s)])2 + Im[S;„(s)]^ 



_ao_ «i / 



(4.16) 



and C'{s) = dC/ds. 

Proceeding in a similar way, this method can be generalized to include non-$ contribu- 
tions and/or the explicit pole structure of the /o(600) and/or additional polynomial terms 
in the Taylor expansion. 



5. Summary 

The procedure we propose for a general un-biased analysis of the e^e" 
sections can be summarized as follows: 



P1P27 cross 



According to Eqs. ( |2.7| )-( pT0 ), the Born (single-photon) cross sections are expressed 
in terms of the three Lorentz-invariant form factors ^1-3. 



The three form factors are decomposed as in (|4.6| ) in terms of a FSR component 
(^fl^) and a two leading SD components {Af^^ and A]^_^l). The FSR component, 
which is present only in the |7r^7r^7) case, is fully determined by the electromagnetic 
pion form factor. 

The vector SD component can be parameterized as in Eq. ( f3.3| ) in terms of the 
effective couplings Fy, gp_^, and g^^,- Here one could use the reference values in 
Table |^ as starting point of the fit, and eventually improve the determination of some 
of these couplings (especially by means of IvrO^o 7) data, where the vector component 
is dominant). 



For the scalar SD component we propose a parameterization of type ( |4.71 ) with a 
resonant part (with or without the /o(600) pole) and a polynomial non-resonant 
term (with one or two free parameters for each channel). In the two |7r7r7) channels 
the parameterization of the non-resonant part can be improved with the inclusion 
of appropriate re-scattering phases, as shown in Eq. ( [4.14 ). These do not involve 
additional free parameters since the Pi in Eq. (|1|) can be determined using the 
precise low-energy constraints on vrvr phase shifts, as outlined in Section 4.1. An 



important consistency check of the whole approach is obtained by the combined fit 
of the two |7r7r7) channels, which should satisfy the isospin decomposition ([1.9|). 



-10- 



Recently, this procedure has been employed in the analysis of KLOE data on e~^e 



7r+7r~-)/ [12 1 with satisfactory results. A more significant test of the method should be 
possible in the near future, with the combined analysis of high-statistics data on both 

e~^e~ — > 7r+7r^7 and e'^e~ — > 7r''7r^7. 
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